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1  INTRODUCTION 


We  consider  only  isotropic  elastic  materials,  and  for  simplicity  thermal  effects  are 
neglected  until  §5.  Under  these  conditions,  the  pressure  is  typically  assumed  to  be 
a  function  of  density  only.  However,  nonlinear  elasticity  theory  predicts  that  the 
pressure  also  depends  on  the  shear  strain,  although  isotropy  implies  this  effect  is 
necessarily  of  second  order;  cf.  Scheidler  [1].  In  §4  we  derive  exact  formulas  for  the 
pressure  and  bulk  modulus  in  a  state  of  uniaxial  strain.  The  effect  of  shear  strain 
can  be  seen  by  comparing  these  results  with  the  corresponding  relations  for  a  state 
of  hydrostatic  stress  (§3).  Our  results  are  based  on  exact  formulas  for  the  speeds 
of  acceleration  waves  (§2).  Applications  to  the  analysis  of  data  from  pressure/shear 
tests  are  discussed  in  §5. 

2  ACCELERATION  WAVE  SPEEDS 

Let  F  denote  the  deformation  gradient  relative  to  the  undeformed  and  unstressed 
state.  The  left  Cauchy-Green  tensor  B  =  FFT  has  principal  values  6,  =  A.2,  where  A,- 
are  the  principal  stretches,  and 


det  F  —  \Jb\  £>2^3  —  —  — ,  p  =  — ,  (2.1) 

P  Po 

where  p  and  p0  denote  the  densities  in  the  deformed  and  undeformed  state.  The 
principal  axes  of  B  are  the  principal  axes  of  strain  in  the  deformed  state.  Since  the 
material  is  isotropic  and  elastic,  these  axes  are  also  the  principal  axes  of  the  Cauchy 
stress  tensor  T,  and  T  is  an  isotropic  function  of  B.  This  implies  that  there  is  a 
single  function  t  such  that  the  principal  stresses  t,-  are  given  by 

=  i(bh  bj,  bk)  =  t(bi,  bk,  bj) ,  (2.2) 

for  any  permutation  i,j,k  of  1,2,3;  cf.  Truesdell  &  Noll  [2,  §48].  It  follows  that  the 
pressure 

P  =  -|trT  =  -|(<i+  <2  +  t3)  (2.3) 

is  a  symmetric  function  of  6j,  62,  b3.  Analogous  results  hold  in  terms  of  the  principal 
stretches  A,-  or  in  terms  of  various  principal  strain  measures,  e.g.,  A<  —  1,  1(6,.  -  1), 
1(1  -  1/6,),  or  in  A,-. 

The  speed  17, •  of  a  longitudinal  acceleration  wave  propagating  along  the  zth  prin¬ 
cipal  axis  of  strain  is  given  by 


nTT 2 

pu<  =2b-dbrXidK 


dti 

din  A,  ‘ 


(2.4) 
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The  speed  Uij  of  a  transverse  (or  shear)  acceleration  wave  propagating  along  the  ith 
principal  axis  of  strain  with  jump  in  acceleration  parallel  to  the  jth.  principal  axis 
(j  /  i)  is  given  by  Ericksen’s  formula: 


if  b^bj, 

b^bj. 


(2.5) 


All  quantities  in  (2.4)  and  (2.5)  are  evaluated  at  the  wave  front.  These  wave  speeds  are 
in  the  deformed  material  (i.e.,  Eulerian);  the  corresponding  Lagrangian  wave  speeds 
axe  obtained  by  dividing  by  A,-.  Proofs  of  (2.4)  and  (2.5)  can  be  found  in  Truesdell 
&  Noll  [2,  §74]  and  Wang  &  Truesdell  [3,  §VI.5].  These  formulae  do  not  require  that 
the  region  ahead  of  the  wave  be  at  rest  or  in  a  homogeneous  state  of  strain.  However, 
when  these  conditions  axe  satisfied,  (2.4)  and  (2.5)  also  apply  to  the  speeds  of  plane 
infinitesimal  sinusoidal  waves;  cf.  Truesdell  &  Noll  [2,  §73]. 


3  HYDROSTATIC  STRESS 

For  a  purely  dilatational  deformation,  we  have 

6,  =  p-2/3  and  t{  =  -p  (i  =  1,2,3).  (3.1) 

In  this  hydrostatic  stress  state,  every  axis  is  a  principal  axis  of  stress  and  strain,  and 
the  pressure  p  is  a  function  ph  of  p  or  p.  Here  and  below,  an  ah”  subscript  denotes 
the  hydrostatic  stress  state.  From  (2.2),  (2.4),  (2.5)x,  and  (3.1),  it  follows  that  for  a 
given  density  p  there  is  a  single  longitudinal  wave  speed  17, •  =  ^  and  a  single  shear 

wave  speed  =  USih ,  and  that 


f  =  uL\~i  tfi;  (3-2) 

cf.  Wang  &  Truesdell  [3,  §VI.5].  A  different  proof  of  this  well-known  result  is  given  by 
Truesdell  &  Noll  [2,  §75].  We  assume  ph  is  a  strictly  increasing  function  of  p.  Then 
(3.2)  implies  the  longitudinal  wave  speed  is  greater  than  the  shear  wave  speed.  With 
the  longitudinal,  shear,  and  bulk  moduli  defined  by 

Lh=pUL^  Gh=pUl  (3.3) 

(3.2)  implies  the  well-known  relation 


Kh  =  Lh-  |  Gh 
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We  use  a  zero  subscript  to  denote  functions  evaluated  at  the  undeformed  and  un¬ 
stressed  state  where  A,-  =  &,•  =  p  —  1;  in  particular,  K0  =  L0-  | G0.  By  (3.4), 


Ph/Ko  «  (p  ~  1)  +  a0(p  ~  l)2 , 


where  the  dimensionless  constant  a0  is  given  by 


_  _L  d2ph 

Q°  2 K0  dp2 


(3.6) 


(3.7) 


4  UNIAXIAL  STRAIN 


For  a  state  of  uniaxial  strain  along  the  1-axis, 

•^i  ==  yfbi =  l//b  A2  =  A3  =  62  =  ^3  =  1  ?  (4-1) 

and  (2.2)  implies  t2  =  t3.  The  principal  stresses  f,-  are  positive  in  tension;  if  <r,-  = 
— f,-  then  <j,-  is  positive  in  compression.  We  use  a  “tt”  subscript  to  denote  uniaxial 
strain  and  consider  only  waves  propagating  along  the  1-axis  into  a  uniaxially  strained 
material.  The  Eulerian  wave  speed  =  Ux  of  a  longitudinal  acceleration  wave  is 
given  by  (2.4)  with  i  =  1,  and  by  (4.1)  we  also  have 


(4.2) 


It  follows  that  a  longitudinal  acceleration  wave  can  propagate  only  if  dax/dp  >  0, 
i.e.,  if  a x  is  a  strictly  increasing  function  of  p,  which  we  now  assume.  By  (4.1),  the 
material  is  strained  iff  p  ^  1  iff  bx  ^  b2.  In  this  case  (2.5)2  and  (4.1)  imply  the  following 
formulas  for  the  Eulerian  speed  U$u  —  U12  of  a  transverse  or  shear  acceleration  wave: 


Gu  =  pUl  =  bx 


tl  —  to  **“  to 


bx  —  b2 
~  f2  _ 
P2-  1  = 


W2 

2  r 


where  r  is  the  shear  stress: 


(4.3) 


r  =  \{ux  -  <r2)  =  |(t2  -  <0  .  (4.4) 

The  Lagrangian  wave  speeds  are  p  Uhu  and  pUSu.  If  p  >  1  the  material  is  in  compres¬ 
sion,  and  (4.3)  implies  that  a  shear  acceleration  wave  can  propagate  only  if  cr1  >  <j2. 
When  p  =  1,  the  results  of  the  previous  section  apply,  and  we  have  Gu|0  =  G0  and 
Ailo  =  L0. 
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From  (4.3)  we  have  the  following  fundamental  formula  for  the  shear  stress  in 
uniaxial  strain: 

r  =  |(^-1  )G„.  (4.5) 

Since  <2  =  t3,  (2.3)  and  (4.4)  imply  the  following  well-known  relation  between  the 
(compressive)  longitudinal  stress  ox,  the  shear  stress  r,  and  the  pressure  pu  in  uniaxial 
strain: 

Pu  =  <*1  -  f  r .  (4.6) 

On  substituting  (4.5)  into  (4.6),  we  obtain  the  following  fundamental  formula  for  pu : 


=  o-i  -  |(p2  -  1)G« 


We  define  the  bulk  modulus  Ku  in  uniaxial  strain  by 


Ku  =  p 


dpu 

dp 


.  dPu  T  dPu 
P  dp  udox  ' 


Then  from  (4.7)  and  (4.2),  we  obtain 


(4.7) 


(4.8) 


Ku  =  Lu- 1 PGU  -  Up2  ~  1  )P^f 

(4.9) 

dG„  1  d 

F„  =  2G„+^  =  --(^). 

(4.10) 

At  p  =  1,  (4.9)  reduces  to  Ku  |0  =  L0  —  | G0  =  K0.  From  (4.9)  and  (4.10)  it  follows 
that  Ku  =  Lu  —  |C7U  for  all  p  iff  Hu  =  0  iff  Gu  —  G0/p2,  but  there  is  no  reason  to 
expect  such  dependence  in  general,  and  thus  no  reason  to  expect  that  Ku  =  Lu  —  | Gu 
except  in  the  limit  of  zero  strain.  Of  course,  by  analogy  with  (3.5)  we  could  have 
defined  Ku  to  be  Lu  —  |C?U,  but  then  (4.8)  would  not  hold.  From  (4.9)  we  see  that 
for  a  state  of  compression,  I<u  <  Lu  -  | Gu  if  Hu  >  0,  and  I<u  >  Lu  -  | Gu  if  Hu  <  0. 
We  assume  that  pu  is  a  strictly  increasing  function  of  p.  Then  any  function  of  p  may 
also  be  regarded  as  a  function  of  ox  or  pu,  and  by  (4.2)  and  (4.8)  we  have 


d_  _  _  _d_  _  _d _ 

P  dp  P  dp  u  do  x  udpu' 


(4.11) 


The  results  up  to  this  point  are  exact.  We  now  consider  some  useful  approximate 
relations.  From  (4.8)  we  have 


Pu/K o  «  (P  ~  1)  +  W  ~  1  )2 4 
where  the  dimensionless  constant  b0  is  given  by 


b0  — 


1  d2pu 


2K0  dp2 


1 

2 


dK,. 


dpu 


(4.12) 


(4.13) 
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For  use  in  (4.12)-(4.13),  note  that  (4.9)  implies 


dKu 

dPu 


Lq  f  dLu 


Kq  \  dal 


8  dG,, 


3  dcr1 


8  Go 
3  Kq 


From  (3.6)-(3.7)  and  (4.12)-(4.13),  we  see  that 

Pu  «  Ph  +  KqCq (p  -  l)2  , 
Pu~Ph  .  Pu  ~  Ph 


Pu  Ph 

where  the  dimensionless  constant  Cq  is  given  by 


Cq  —  bo  —  o,q  — 


1  (dKu 

2  V  dPu 


Co(p~l), 


dKh 


dPh 


(4.14) 

(4.15) 

(4.16) 

(4.17) 


On  comparing  (4.16)  with  equation  (4.6)  in  Scheidler  [1],  we  find  that  Cq  is  also  given 
by 


Co=3 


2  (dGh 
dPh 


2l 
Ko t 


(4.18) 


5  DISCUSSION 

The  longitudinal  stress  ^  as  a  function  of  p  in  uniaxial  strain  can  be  obtained 
from  normal  plate  impact  tests.  Then  the  relation  (4.6)  (which  does  not  rely  on  the 
assumption  that  the  response  is  elastic)  is  typically  used  to  determine  the  pressure 
pu  in  uniaxial  strain  given  some  assumptions  on  the  shear  stress  r,  or  to  determine  r 
given  some  assumptions  on  pu.  It  is  often  assumed  that  pu{p)  is  equal  to  the  pressure 
Ph{P)  in  a  state  of  hydrostatic  stress  at  density  p  (or  to  some  appropriate  modification 
of  Pf,  to  include  thermal  effects  in  the  shocked  state).  Such  an  approximation  neglects 
the  effects  of  shear  strain  (or  shear  stress)  on  pu.  That  this  effect  may  be  significant 
in  ceramics,  geologic  materials,  and  polymers  has  been  emphasized  by  Gupta  [4] 
and  Conner  [5].  These  materials  can  sustain  relatively  large  elastic  shear  strains 
(compared  to  metals),  although  for  polymers  viscoelastic  effects  should  also  be  taken 
into  account.  Only  elastic  response  is  considered  here.  Then  (4.15)  implies  that  pu(p) 
differs  from  Ph(p)  by  a  term  of  order  (p  —  l)2  unless  Cq  =  0,  which  is  generally  not 
the  case.  If  Cq  and  ph{p)  are  known,  then  (4.15)  provides  an  approximation  to  pu  to 
within  an  error  of  order  (p  —  l)3.  The  relative  error  in  approximating  pu  by  p^  is  of 
order  p  —  1  and  can  be  estimated  by  using  (4.16). 

In  a  pressure/shear  (or  oblique  plate  impact)  test,  a  longitudinal  wave  propagates 
into  the  undeformed  material,  bringing  it  to  a  state  of  uniaxial  strain,  and  a  slower 
shear  wave  propagates  into  this  uniaxially  strained  material.  These  tests  yield  both 
cr-i(p)  and  the  shear  wave  speed  Us<u  (and  hence  G„)  as  a  function  of  p  or  ax.  If 
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the  shear  wave  travels  at  the  acceleration  wave  speed,  then  (4.5),  (4.7),  and  (4.9) 
provide  exact  formulas  for  the  shear  stress  r,  the  pressure  pu,  and  the  bulk  modulus 
Ku  in  uniaxial  strain  as  a  function  of  p  or  <jx .  These  formulas  appear  to  have  gone 
unnoticed,  however.  Instead,  it  is  usually  assumed  that  Ku  «  Lu  —  \GU.  This 
approximate  relation,  together  with  (4.8),  is  then  integrated  to  give  pu  as  a  function 
of  p;  cf.  Gupta  [4,  6],  Conner  [5],  and  Aidun  &  Gupta  [7].  For  fused  silica  in  the  strain 
range  0  <  p  —  1  <  0.076,  the  response  is  elastic  and  the  shear  wave  speed  decreases 
with  />;  cf.  Conner  [5].  In  this  strain  range  the  shear  wave  is  an  acceleration  wave  (cf. 
also  Abou-Sayed  &  Clifton  [8]),  so  we  may  apply  the  results  of  §4.  Using  (4.9)  and 
Conner’s  data,  we  find  that  at  a  strain  of  p  —  1  =  0.076  the  estimate  Ku  L.  - 1 G. 
is  low  by  about  29%. 

Whether  the  shear  wave  in  a  pressure/shear  test  is  an  acceleration  wave  or  a 
shock  wave  depends  on  the  nonlinear  elastic  response  of  the  material;  cf.  Davison 
[9].  The  shear  modulus  Gu  in  §4  is  defined  in  terms  of  the  acceleration  wave  speed 
or  equivalently,  in  terms  of  the  speed  of  a  plane  infinitesimal  sinusoidal  shear 
wave;  cf.  §2.  If  a  shear  shock  with  speed  U  can  propagate  in  the  uniaxially  strained 
material  and  if  we  set  G  =  pU 2,  then  the  formulas  in  §4  hold  approximately  when 
Gu  is  replaced  with  G.  Also  note  that  if  U  >  USu  (as  standard  stability  arguments 
would  imply),  then  G  >  Gu,  and  (4.5)  and  (4.7)  imply  that  r  <  |(/52  —  1)(5  and 
pu  >  <Ti  -  |(/52  —  1  )G  in  compression  ( p  >  1). 

We  conclude  with  a  brief  discussion  of  thermodynamic  effects,  which  have  been 
neglected  up  to  this  point.  If  a  thermoelastic  material  conducts  heat  by  Fourier’s  law 
[respectively,  is  a  nonconductor],  then  a  longitudinal  acceleration  wave  propagates 
at  the  isothermal  [respectively,  adiabatic]  wave  speed.  However,  the  formula  (4.3) 
for  the  speed  of  a  shear  acceleration  continues  to  hold  in  either  case;  cf.  Bowen 
&  Wang  [10].  In  fact,  it  can  be  shown  that  (4.3)  holds  even  if  heat  conduction 
is  governed  by  Cattaneo’s  equation,  which  prohibits  instantaneous  propagation  of 
thermal  disturbances.  Thus  the  formulas  (4.5)  and  (4.7)  for  the  shear  stress  and  the 
pressure  continue  to  hold.  In  particular,  they  are  valid  when  the  state  of  uniaxial 
strain  has  been  achieved  by  shock  loading. 
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